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Abstract 



We study the order parameter distribution P(q) in the 4d Ising 
spin glass with ±J couplings in a magnetic field. We also compare 
these results with simulations for the infinite ranged model (i.e. SK 
model.) Then we analyse our numerical results in the framework of 
the droplet picture as well as in the mean field approach. 
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This work is devoted to the study of spin glasses in presence of a magnetic 
field. During the last ten years, a large number of works has been devoted 
to the study of spin glasses [Q]. One of the problems which still remains 
unsolved is to understand the effect of a magnetic field on the spin glass 
phase. The mean-field theory predicts that the spin glass phase will survive 
to the application of a magnetic field below the de Almeida-Thouless critical 
line (AT line) ||. In the most general case, the main effect should be the 
destruction of a large number of equilibrium states with a reshuffling of the 
free energies for the remaining ones. To our knowledge, even at the mean field 
level, a numerical test of the theoretical predictions of the replica symmetry 
breaking solution with magnetic field has never been done. Such a test would 
be interesting because it would give support to the Parisi ansatz as a correct 
solution to mean- field theory [||]. 

For the short-range models case, there is still much controversy. In fact, 
there is no precise theoretical prediction. A usual e expansion near dimension 
6 will run in trouble because it is not known how to find a non trivial fixed 
point [|J. Phenomenological models like those developed by D. S. Fisher and 
D. A. Huse [§] predict that the spin glass phase disappears for a finite mag- 
netic field. But this result is a consequence depending on some assumption 
on the real nature of the low temperature spin glass phase. The most recent 
studies of the AT line were done in Monte Carlo simulations. In these works, 
the main points of interest were the curves of constant non-linear suscepti- 
bility in the h — T plane || or its divergence when approaching the AT line 
using finite-size scaling methods 0. Nevertheless the first approach is very 
indirect and the second one could be plagued by strong corrections to the 
simple scaling. 

In this work we have tried to understand the spin glass with magnetic 
field by studying the P{q) order parameter function obtained by means of 
the Monte Carlo method using the heat bath algorithm. We will first present 
a brief discussion on the general theoretical predictions of this problem and 
numerical results for the mean-field theory. Then we will present the numer- 
ical results obtained for the id ± J Ising spin glass. It is well established that 
this model has a finite T c and it has been intensively studied || . 

The (i-dimensional Ising spin glass model of interest, with ±J couplings, 
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is defined by the following hamiltonian 

H = - ^2 JijViVj - h^ai (1) 

(m) ' 

The couplings are quenched variables with equal distribution values ±1. 
The interaction is restricted to nearest neighbors and h is the magnetic field. 
The Ising spins take two possible values ±1 and live in a <i-dimensional 
hypercubic lattice with periodic boundary conditions. It is very useful to 
consider discrete couplings in the hamiltonian because this speeds up the 
updating of the spins in the Monte Carlo numerical simulation and their 
discreteness should not be relevant for the physics at least for not too low 
temperatures. In the limit d — > oo, one expects to converge to mean-field 
theory, i.e. the SK model []1U|] . In the SK model, all spins interact among 
them and the couplings are normalized by a factor 1/vJV where N is the 
number of spins. 

We consider two identical copies of the system eq.(|lj), i.e with the same 
realization of the bond disorder [|ll[]. Let us call them {<jj} and \ri\. The 
overlap Q among the two copies is defined by 

Q = I E '<« ( 2 ) 

i 

from which we can construct the order parameter function P(q) 

P(q) = (S(q-Q)) (3) 

where (...) and (...) mean the usual statistical Gibbs average over configura- 
tions and the average over the quenched disorder respectively. 

In mean-field theory, below the critical temperature and at zero magnetic 
field, there exist an infinity of equilibrium states, all of them having a different 
statistical weight. Furthermore, all these states have zero magnetization and 
no state is particularly selected if we apply a magnetic field. In fact, an 
infinity of equilibrium states still remain and the spin glass phase survives 
in a magnetic field. At zero magnetic field the order parameter distribution 
P(q) is symmetric under the exchange q — > —q. The magnetic field breaks 
this symmetry and the P(q) is expected to be non-zero only if q > 0. This 
means that half of the states have been suppressed by the magnetic field. 
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Close to T c , the free energy of the SK model with magnetic field can be 
approximate by 

f = rJ2Ql-l Tr Q 3 - y 2 E QU - h2 J2 Q* ( 4 ) 

a<b a<b a<b 

where {Q a b] 1 < a, 6 > n} is the order parameter and n the number of repli- 
cas. The equilibrium solution of the free energy in the limit of infinite order 
of replica symmetry breaking gives a function q(x) defined in the interval 
(0, 1) ||. This function q(x) is the analytical continuation of the matrix 
in the limit n —>■ 0. Close to T c one gets a q(x) with two plateaus in the 
regions < x < x min and x max < x < 1 (x min < x max ) with respective 
values q min and q max (q min < q ma x)- Between x min and x max , q(x) increases 
with x. q max is nearly independent on the field but q min increases with a 
power of h smaller than 1 (q m i n ~ h*). Using the static chaos approach to 
spin glasses it has been suggested that the effect of the magnetic field is 



the progressive suppression of all equilibrium states a such that their overlap 
< Qmin V/3. This corresponds to cutting some branches of the ultrametric 
tree. Using the relation P(q) = one finds that P{q) is given by a 

continuous part Po(q) which is non zero inside the interval (q m in,qmax) and 
two singularities at the extremes of this interval 

P(q) = P (q) + a5(q- q min ) +bS(q- q max ) (5) 

We have simulated the SK model at T = 0.5 and h = 0.3 (the corre- 
sponding field at the AT line at that temperature is h ~ 0.57). We simulated 
three different sizes N = 320, 1048 and 3200. For these sizes we were able 
to reach equilibrium for near all samples after 100000 Monte Carlo steps for 
the largest size. Then statistics was collected over several hundred thou- 
sands of Monte Carlo steps. The main source of fluctuations comes from 
the finite number of samples because the P{q) is strongly non self-averaging. 
Self-averageness is restored when the AT line is reached by increasing the 
field. The number of samples is 500, 30 and 20 respectively. Even though the 
numbers of samples are small, they are large enough to show the qualitative 
behavior of the P(q)- The results are shown in Fig. 1. The P{q) begins to 
display two singularities for sizes of several thousands spins. For smaller sizes, 
we only found one peak plus a long tail which extends down to the region 
of negative overlaps. According to the Parisi solution to mean-field theory, 
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q m i n should match the correct value at infinite order of replica symmetry 
breaking in the infinite-size limit. We can then compute the position of both 
singularities, at least at first order of replica symmetry breaking and this 
gives q min = 0.45, q max = 0.63. This is in agreement with our numerical re- 
sults. Other good estimates for q min and q max are also obtained using the PaT 
(Parisi-Toulouse) hypothesis JT3|] (which is a very good approximation at least 
close to T c ). This approximation predicts that q min (h,T) = q(h, T^rih)) and 
q m ax{h,T) = q(IiAT(T),T) where TatW and }iat(T) are the equations for 
the AT line. Computing these values, one gets q max ~ 0.64 and q m i n ~ 0.437. 
This is also in agreement with our simulations. 

Now we return to the 4d case. There exists two possible scenarios that 
we want to compare. First, from the droplet models |5j it is expected that all 
excitations of droplet of sizes larger than a certain length £ will be suppressed 
by the field. The dependence of this correlation length in function of the 
magnetic field is given by 

i~{q E Ah 2 )^ (6) 

with qEA the Edwards- Anderson order parameter and 9 the thermal exponent 
which gives the characteristic energy scale L e of droplet excitations of typical 
size L. This exponent 9 should be approximately (as emerges from 
numerical studies of chaos in spin glasses 0, [L5|.) For sizes much larger 
than £ it is expected that the P(q) will be strongly peaked around a unique 
value of q. 

We have simulated L = 3, 5, 6, 8 in a 4d lattice with periodic boundary 
conditions with ±J couplings. Simulations were performed at T = 1.2 (~ 
0.6 T c ) and h = 0.4. The number of sample are 320, 128, 100, 50 respectively. 
From finite-size scaling studies 0, we expect to be within the spin glass 
phase if there is an AT line. It is not very difficult to reach the equilibrium 
in case of L = 3, 5. 100000 Monte Carlo steps were enough after a slow 
cooling procedure. Now we will try to convince the reader that we effectively 
thermalized for L = 6,8. To this end we performed a simulated annealing 
of half a million of Monte Carlo steps from the high to the low temperature 
phase at constant magnetic field. After that, statistics was collected over 
the next half a million Msteps. During this collecting, we computed the four 
moments of the P{q) distribution which show no apparent drift in time. To 
increase the statistics we simulated in parallel eight identical copies of the 
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system computing the four overlaps among four different pairs at each Monte 
Carlo step. The fact that the magnetic field tries to align the spins helps 
in the thermalization procedure. This is the reason why we were able to 
thermalize over a scale of time of half a million of Monte Carlo steps which 
would be probably insufficient at zero magnetic field. Figure 2 shows the 
numerical results for the P(q). We can immediately notice that there is no 
singularity at q — q min if we compare to the previous figure for the SK model. 

Looking at this results it is difficult to draw a definite conclusion on what 
is the correct scenario in 4d Ising spin glasses. Two facts are interesting to 
point out. The first one is the existence of a long tail for sizes up to L = 6 
which extends down to negative overlaps. So, P(q = 0) is finite which means 
that reversal of compact domains of characteristic size L are still present with 
a finite probability. Within the droplet model, we can estimate how domain 
excitations of typical size L are suppressed by the magnetic field. The effect 
of the magnetic field depends on the regime in which the system is, either 
L>>£orL<<£, £ being given by eq.(|). We can estimate the value 
of £ by using numerical simulations of static chaos A typical value of 
order 5 is obtained. If L >> £ we expect that the droplets excitations of size 
L are suppressed with a factor ~ exp(— j3x(L) d h 2 ) respectively to the case 
h — (x being the linear susceptibility.) So, in this regime, the tails would 
be suppressed. Unfortunately, we are in the regime where £ ~ L. As a lower 
bound, when L « £, tails are suppressed with a factor ~ exp(—(L)2h 2 ). 
In our range of sizes, this factor is of order 10" 1 which is smaller than what 
we can see on our plots (P(q = 0) being 0.3 at zero magnetic field, we would 
expect P(q = 0) ~ 0.03 [§.) The second fact regards the absence of a second 
peak of P(q) at q m in- Presumably, such a peak could appear for larger sizes. 
In the case of SK model, the q m i n peaks already arise for size of order 1000 
spins, as oppose to the 4d case. One possible reason for such a difference 
reside in the fact that, for the 4d case, we can be very close to the AT line. 
A second reason is that we can surely expect stronger finite size effects than 
in the mean-field case. For instance, the singularity in P(q) for q = q max is 
less pronounced as can be seen in numerical simulation |5J. 

In order to reach more definite conclusion, we need to study larger sizes 
lattices. In practice, such a task is very difficult because for larger sizes 
lattices we are not able to thermalize. In fact, we have performed numerical 
simulations for L = 10 and L = 12. Despite that these are non equilibrium 
results, interesting hints can be obtained. In Figure 3 we show the P(q) 
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distribution. Starting from uncorrected configurations, the overlap among 
two copies grows with time. In several cases it remains stacked in a value of q 
close to 0.4 giving two singularities for the P(q) distribution. This indicates 
that we are in the good region in the h — T plane in order to test if there 
exists a spin glass phase. 

Still, in order to have a more definite conclusion, we need to take advan- 
tage of new numerical simulation techniques like the simulated tempering 



lq| . This method has revealed much effective for the 2c? |17j and 3d |18 
Ising spin glasses. A work using such techniques is under progress. 

Summarizing, we have studied the 4d Ising spin glass with magnetic field. 
For comparison, we have also simulated the SK model. This is also a test 
of the Parisi solution to mean-field theory and our numerical results are 
in agreement with it. In the 4d case we present results for lattice size up 
to L = 8. Then we tried to interpret them in the mean- field picture and 
the droplet one. It seems that the effect of the magnetic field is weaker than 
what droplet picture predicts. Non thermalized results for larger sizes suggest 
that fully equilibrated simulations should be able to select in a definite way 
between these two pictures. We hope that using numerical techniques like 
simulated tempering should be able to decide the question in the near future. 



Acknowledgments 

We thank D. Huse and G. Parisi for very interesting suggestions. F. R. has 
been supported by the EC fellowship B/SC1*/915198 and the INFN. M. P. 
has been partially supported by the EU Science grant SC1*0394 and EU 
Capital Humain et Mobilite grant - Euronetwork ERB 4050PL92 0982/0035. 



S 



References 

[1] M. Mezard, G. Parisi and M. A. Virasoro, "Spin Glass Theory and 
Beyond" World Scientific, (1987); K. Binder and A. P. Young, Rev. 
Mod. Phys. 58 (1986) 801; K. H. Fischer and J. A. Hertz, "Spin Glasses" 
Cambridge University Press, (1991). 

[2] J. R. de Almeida and D. J. Thouless, J. Phys. All (1978) 983. 

[3] G. Parisi, J. Phys. A 13 (1980) 1101; J. Phys. A 13 (1980) 1887. 

[4] A. J. Bray and S. A. Roberts, J. Phys. C 13 (1980) 5405. 

[5] D. S. Fisher and D. A. Huse, Phys. Rev. Lett. 56 (1986) 1601; Phys. 
Rev. B 38 (1988) 386. 

[6] E. R. Grannan and R. E. Hetzel, Phys. Rev. Lett. 67 (1992) 907. 

[7] D. Badoni, J. C. Ciria, G. Parisi, J. Pech, F. Ritort and J. J. Ruiz, 
Europhys. Lett. 21 (1993) 495; J. C. Ciria, G. Parisi, F. Ritort and J. 
J. Ruiz, J. Phys. I. (France) 3 (1993) 2207. 

[8] J. D. Reger, R. N. Bhatt and A. P. Young, Phys. Rev. Lett. 16 (1990) 
1859; G. Parisi and F. Ritort, J. Phys. A 17 (1993) 3179; J. C. Ciria, 
G. Parisi and F. Ritort, J. Phys. A 17 (1993) 3179. 

[9] E. Marinari, G. Parisi and F. Ritort, preprint |cond-mat 931004T . 

[10] D. Sherrington and S. Kirkpatrick, Phys. Rev. B 17 (1978) 4384. 

[11] A. P. Young, Phys. Rev. Lett. 51 (1983) 1206. 

[12] C. De Dominicis and A. P. Young, J. Phys. A 16 (1983) 2063. 

[13] G. Parisi and G. Toulouse, J. Phys. Lett. (France) 41 (1980) L361. 

[14] I. Kondor J. Phys. A 22 (1989) LI 63. 

[15] F. Ritort, submitted to Phys. Rev. B. 



9 



[16] B. A. Berg and T. Celik, Phys. Rev. Lett. 69 (1992) 2292; B. A. Berg 
and T. Neuhaus, Phys. Lett. B 67 (1992) 249; E. Marinari and G. Parisi, 
Europhys. Lett. 19 (1992) 451. 

[17] W. Kerlerm and P. Rehberg, preprint |cond-mat 9402049 . 

[18] E. Marinari, G. Parisi and F. Ritort, in preparation. 



10 



Figure Captions 



Fig. 1 P(q) for the SK model at T = 0.5, h = 0.3. The error bars are of 
order 20% for N = 3200 and 15% for N = 1408 and less than 5% for 
N = 320. The symbols are a guide to the eye. 

Fig. 2 P(q) for the Ad Ising spin glass. Error bars are smaller than 15% in all 
cases. The symbols are a guide to the eye. 

Fig. 3 Non thermalized P(q) for the Ad Ising spin glass. 
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